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A theoretical analysis for the dynamic performance characteristics of 
laminar flow, capillar/orifice compensated hydrostatic journal bearings is 
presented. The analysis considers in detail the effect of fluid compressibility 
in the bearing recesses. At high frequency excitations beyond a break frequency, 
the bearing hydrostatic stiffness increases sharply and it is accompanied by a 
rapid decrease in direct damping. Also, the potential of pneumatic hammer 
instability (negative damping) at low frequencies is likely to occur in 
hydrostatic bearing applications handling highly compressible fluids. Useful 
design criteria to avoid undesirable dynamic operating conditions at low and 
high frequencies are determined. The effect of fluid recess compressibility is 
brought into perspective, and found to be of utmost importance on the entire 
frequency spectrum response and stability characteristics of hydrostatic/hybrid 
journal bearings. 


INTRODUCTION 

The dynamic performance of hydrostatic and hybrid journal bearings (HJBs) 
has been a subject of growing theoretical interest. In the past, conventional 
applications of HJBs were restricted to support heavy loads at zero or low 
speeds and also to low operating eccentricities required for precise machine 
tools. However, industry's expanding demands for higher speed applications and 
the ability of bearings to support dynamic loads have brought the need to 
consider and study in detail the dynamic performance of these type of bearings. 

Davies and Leonard (1969) first studied the dynamic behavior of HJBs, and 
derived under restricted conditions closed form solutions for sjnall motions 
about the concentric position. The limited stability characteristics of hybrid 
bearings have been confirmed experimentally by Davies and Leonard (1970, 1971), 
and Leonard and Rowe(1973). M.K. Ghosh (1979) and V.N. Singh(1979) have also 
presented numerical analyses for the static and dynamic performance 
characteristics of HJBs. For best dynamic performance, the analyses show that 
optimum direct damping is attained at low recess pressure ratios and reduced 
recess to bearing area ratios. Based on the concentric thin land theory, Rowe 
(1980) presented coefficients for hydrostatic and hydrodynamic stiffness and 
squeeze damping for capillary, orifice and constant flow control HJBs. The 
mechanism of whirl instability for hybrid operation is explained, and 
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2 2 
n d /4. Orifice area [m ] 

Recess circumferential length [m]. 

radial clearance [m] 

damping coefficients [Ns/m] 

Orifice discharge coefficient 

2 R . Bearing diameter [m]. 

Orifice and capillary diameter [m] 

oo/o^. dynamic frequency ratio. 

H/c = 1+Ae^cosG + As^sinG Dimensionless film thickness. 

cosG, sinG 

Recess depth [m] . 

Stiffness coefficients [N/m] 

Bearing axial width, recess axial length [m]. 

Number of recesses on bearing. 

2 

Fluid pressure [N/m ] . 

2 

Ambient and external supply pressures [N/m ] 

(P -P )/(P -P ). recess pressure ratio. 

LG. So. 

3 

Q y/c (P -P ) Dimensionless recess flow rate into film 
r 3 a lands 

■]2(P s -P r )/ p. Fluid speed through orifice [m/s]. 

'll/op) . Orifice Fluid speed at sonic conditions [m/s] 

’lA /( |3 8nyl c ) . Capillar fluid speed at sonic conditions 

(H +H)lb+V Total recess volume, volume of orifice supply 
line [ in ]. 

Inertial coordinate system, 
restrictor control parameter. 

reduced damping factor due to compressibility, eqn. (14) 

2 

Fluid compressibility factor [m /N], 
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Ao p/2/c // p(Ps-Pa) Design Orifice parameter. 

2 3 

A / ( 8 Tl 1 c ). Design capillary parameter. 

2 

yaA^/c (Ps-Pa) Recess frequency parameter 

V 0(Ps-Pa)/clb. Compressibility parameter 
2 2 

yS2l /c (Ps-Pa) Rotational speed parameter. 

2 2 

yoit /c (Ps-Pa) Frequency parameter at film lands. 

Excitation or whirling frequency [1/sec] 

Break frequency [1/s], eqn. (15). 

Rotational speed of journal [1/s] 

K / 2C . Whirl frequency ratio. 
xy xx 'l j 

Dimensionless perturbed eccentricities 
2JI/N; b/R 
50-5$; (L-l)/D 

cot Dimensionless time coordinate 

(x,y)/R. Dimesionless circumferential and axial coordinate 
2AEL 

Circumferential flow factor 

50 A0 D 


Subscripts: 
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a 


Superscripts: 
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refers to zeroth order solution 
or no compressibility solution. 

refers to bearing recesses 

refers to first order solution 

refers to n-th recess 
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it is shown that the whirl frequency ratio is equal to 0.5. Nonetheless, the 
introduction of hydrostatic stiffness into the rotor-bearing system raises the 
frequency of whirl onset. 

Rohde and Ezzat (1976) presented an analysis which accounts for the 
lubricant compressibility in the bearing recesses and supply line. The results, 
also supported by the analysis of Ghosh and Viswanath (1987) and Ghosh et al. 
(1989), show that when the journal is subjected to high frequency excitations 
the bearing stiffness and damping characteristics can change drastically. The 
behavior is determined by a "break" frequency beyond which the bearing direct 
stiffness increases sharply and it is accompanied by a rapid decrease in direct 
damping. The analyses cited are based on full numerical solutions to the 
Reynolds equations, and thus, their results can not be easily used in practice. 

The present work introduces a theoretical analysis for the dynamic 
performance characteristics of capillar/orifice compensated centered hydrostatic 
journal bearings. The analysis considers in detail the effect of fluid recess 
compressibility. Useful design criteria to avoid undesirable dynamic operating 
conditions at low and high frequencies are determined. Conditions to avoid 
complete loss of damping at small frequencies are also presented. 

The relevance of the analysis is immediately recognized since the performance 
characteristics of HJBs are relatively independent of the static journal center 
displacement for small to moderate values of the eccentricity vector. 


ANALYSIS 


The geometry of a hydrostatic bearing and reference coordinate systems 
are shown in Figure 1. A lubricant at external pressure P is supplied through 
either a capillary or orifice to each recess, flows into fhe film lands, and 
then exits the bearing. The performance characteristics of HJBs are governed by 
momentum and continuity considerations at the film lands as well as flow 
continuity through the recesses. The former define the pressure and flow fields 
over the thin film region, whereas the latter provides values for the pressure 
in the recesses. 


The continuity equation at each recess determines the global flow 
balance between the flow through the orifice (Q ) or capillary (Q ) restrictors, 
the recess outflow into the film lands (Q ), anc? the time change of fluid mass 
within the recess volume. The recess continuity equation is then given as 


9V 9P 

Q. = Q + — + V 0 — 
i r r 


9t 


9t 


( 1 ) 


where 0 represents the fluid compressibility factor. The inlet flow (Q^) through 
the external restrictor depends on the nature of the control device, i.e. 


for orifice: 


Qi = Qo = Cd 


nd 2 
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(2a) 
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for capillar: 


nd 


Q i " Q c 


128yl 


(P - P ) 
v s r' 


(2b) 


For small amplitude motions of frequency w about the concentric position, the 
land film thickness is described by the real part of the expression: 


(3) 


h = H/c = 1 + e 1T ( Ae h + Ae h ) 

xx y y 

h x = cos 0 , h y = sin 0 

T=wt, Ae , Ae « 1 
x y 

The film land and recess pressures are expressed as the sum of a zeroth and 
first order complex fields describing the equilibrium centered position and the 
perturbed dynamic motion, respectively. Hence, in dimensionless form 


i x 

p = p + e (Ae p +Ae p) 
r *o x *x y *y' 

i T 

p = p+ e (Aep+Aep) 
*r *ro v x *rx y v ry J 


(4) 


The linear perturbation method used determines the recess continuity equation to 
be accordingly divided into zeroth and first order expressions as 

zeroth order recess continuity equation: 


a) capillary: 

b) orifice: 


5 c 


p ro^ q ro 


5 i ( 1— p ) = q 
o v ‘ro 7 M ro 


(5a) 

(5b) 


and the dimensionless film flow into the lands is given by (Davies, 1969, Rowe, 
1980) as 

59 


q 


ro 


1215 Pro 


( 6 ) 


5 and 5 are the capillary and orifice design parameters that help to 
determine the restrictor geometry for the desired operating pressure ratio. 

first order recess continuity equation at nth recess: 


a)capillary: 

-Pra" ' 5 c * 1 h *2 1 - Ira 

(7a) 

b)orifice: 

5 2 

n r o • \ \ i n 

‘Pra 1 — + 1 h h J ' "ra 
/q ro 

(7b) 


a = x,y; n=l,2, . . . . ,N 
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and q n is the perturbed land flow on the closure surrounding the nth recess. 
For laminar flow, the film flow is equal to (Rowe, 1980, SanAndres, 1989a) 


p ro 

9 i 

_ n+1 

h d0 + 
a 

0 

n 

2 S© 

n 

. . . . hydrostatic 

axial 

flow 

2 A£ v 

P rw * 

12 A£ ra 

L/D 

+ 

[ 2 p n 

n+1 n-1 . 

- p - p ] 

. . . . hydrostatic 

circ. 

flow 


( 8 ) 

hydrodynamic flow 
squeeze film flow 

0 = nS9 
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Collecting similar terms, equations (7) are reduced to 
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where z is a control parameter, first introduced by Rowe (1980), and given as 
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for orifice and capillary restrictors, respectively. Analytical solution to the 
system of equations (9) is easily accomplished(see Davies and Leonard, 1969). 
The resulting dynamic recess pressures are equal to 
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where f=oo/c^ is a frequency ratio, c*^ is named a "break" frequency, and its 
relevance on the bearing dynamic force response will be discussed latter. The 
break frequency is given by the expression: 


% = 


TtDc 3 

(z+l+2ysin ( it/N) ) 

3N(L-l)p|SV r 


( 11 ) 


Dynamic force coefficients are calculated by integration of the first order 
pressure field over the journal surface, 


K «e + lM c m3 
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a, (3 = x, y 


and it is recalled that h =cos 0 and h =sin 9, respectively. After some 
algebraic manipulations, it is found tKat the stiffness and damping 
coefficients are given by 
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where K xxq , C , K q are t ^ e dynamic force coefficients in the absence of the 
fluid compressibility effect. These coefficients are identical to those 
presented by Rowe (1980) and presented in Table 1. 

The parameter a is denoted as a reduced damping factor due to 
compressibility at zero frequency and given by 


3 N V r (P r -P a ) (3 

c = 2 

JtLDc [z+1+2 ysin t n/N) ] 


(14) 


The dynamic force coefficients given in equations (13) depend strongly 
on the frequency ratio f= oVoig. Dynamic operation at f=l results in a reduction 
of 50% in direct damping and a considerable increase in hydrostatic stiffness. 
The "break" frequency equation (12), is rewritten in terms of the restrictor 
characteristics as 


a) orifice HJB: 


“b = 


A v 
o so 


V v z 
r o o 


(z o +l+2Ysin 2 (n/N)) 


(15a) 


b) capillary HJB: 

A v „ 

= c (z +l+2Ysin 2 (jt/N)) (15b) 

V z c 

r c 

where v is the fluid velocity at sonic conditions across the orifice or 
capillar, and v is the fluid speed through the orifice restrictor. Appropriate 
design of HJBs requires the bearing to operate under dynamic conditions well 
below the break frequency ocl. As a rule of thumb, it is suggested that operation 
above f=0.01 should be avoided. High break frequencies can only be achieved on 
HJBs with shallow recesses (small overall recess plus supply line volumes) and 
large diameters for the orifice and capillars. For orifice compensated HJBs, a 
large ratio of sonic speed to fluid orifice speed is desired; while for capillar 
HJBs, the restrictor length should be kept to a minimum. Of course a negligible 
value of the fluid compressibility factor B is sufficient to insure adequate 
dynamic operation without stiffness hardening and reduced damping. It is 
mentioned that fluids of large viscosity and HJBs with small axial recess 
lengths are more likely to present low values of the break frequency o^. 

Dynamic operation at very high frequencies determines a complete loss of 
damping while the hydrostatic stiffness attains its limiting value. As <o ®, 
from equation (13a) we obtain 
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Note that the asymptotic value of direct stiffness depends solely on the 
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compressibility factor (3 and the recess volume, and it is independent of the 
fluid properties or external supply conditions. 

Although dynamic operation at frequencies well below is highly 
desired, the possibility of "pneumatic hammer instability" in a HJB should be 
carefully analyzed. This phenomenon characteristic of highly compressible fluids 
such as in gas-lubricated hydrostatic bearings (Lund, 1965) occurs at low 
frequencies and determines negative values of the direct damping coefficients. 

At zero frequency, from equation (12b), the value of direct damping is given as 
C xx = C xxo )• Thus, null or negative damping is avoided if the reduced 

damping ractor £s small, a,«l, or from equation (14), if 

N V (P -P ) 3 1 

— - « - [z+1+2y sin Z (n/N)] (17) 

It LDc 3 

which requires a low ratio of recess volume to film land volume, and small 
products of recess pressure and compressibility factor 3* This condition possesses 
a severe restriction on the recess volume and the maximum allowed recess 
pressure for highly compressible fluids such as liquid hydrogen. Redecliff and 
Vohr (1969) have recommended that the total volume in the feeding recesses be 
less than 5% the total volume in the bearing film in order to avoid pneumatic 
hammer effects. The relation derived in the present analysis gives a better 
though more stringent criteria. Equation (17) requires the recess volume to be 
low, and thus, the recess depth must be kept to a minimum consistent with 
providing uniformity of pressure over the recess. 

For hybrid operation, the whirl frequency ratio is a fundamental 
property of the bearing( independent of the rotor system) which determines the 
maximum shaft operating speed before whirl instability occurs. The whirl ratio <f> 
is defined as 

K 1 

♦ (18) 

52 C 2 (1-a ) 

Thus, in a centered HJB the whirl ratio <f> is at most equal to 0.50, and it 
is substantially increased if the reduced damping factor a is large. From a 
rotordynamics point of view, hybrid bearings offer no advantage in relation to 
plain journal bearings. Additionally, if fluid recess compressibility is large, 
it causes the limited applicability of hybrid bearings to be further restricted. 

The conditions stated are of considerable importance in RJBs used for 
cryogenic applications such as those handling liquid hydrogen. For LH2, the 
compressibility factor 3 is rather large while its viscosity is very small. 
Typical bearing supply pressures exceed 6.9MPa(1000 psig) in order to attain the 
required levels of load capacity. In this case, even if recesses of low depth 
are designed, it is very likely that the reduction in damping due to the fluid 
compressibility is large. 


EXAMPLES 

The following sample case serves to show the effect of recess 
compressibility on the dynamic force coefficients of a typical capillary HJB. 
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The example has been taken from an application presented by Rhode(1976). The 
geometric and operating parameters for the bearing are as follows; 

N=8, L=35.56 mm(1.4in), D=71.12 mm(2.8in), 

1=30.48 mm(1.2in), b=22.84 mm(0.9in), 

c=0.0254mm(0.001in) , Hr/c=10,30,50, 

d =0. 325mm, 1 =19. 1mm, 
c c 

y=0.0345 Ns/m 2 , 3=5.8 10~ 10 m 2 /N, 

P -P = 689.5 KPa(1000 psig), p =0.3 
8 S = a 1000 rpm (16.66 Hz), ro 

For the sample bearing, the static direct stiffness and damping _ 

coefficients at zerg frequency and no compressibility are equal to K =5.20 10 / 
N/m and C xxq = 2.5 10 Ns/m, respectively. For 3 different values of tiie°recess 

depth, Figures 2 and 3 show the variation of K and C for increasing values 

of the frequency ratio f. The calculated valuesof the X ^reak frequency and 
reduced damping factor are given by 

Hr/c=10; X,=0.0044, 0^=2281 Hz, a =0.001448, 

Hr/c=30; XT=0. 0124, av= 809 Hz, a C =0. 004081, 

Hr/c=50; >^=0.0204, <*£= 492 Hz, o£=0. 006710, 

Figure 2 shows that the direct stiffness increases dramatically as the 
excitation frequency grows. Note that the greater the break frequency(smaller 
recess depth), the greater is the asymptotic value of the stiffness coefficient. 
Figure 3 shows the damping coefficient to decrease as the excitation frequency 
increases. Note that the damping curves are almost identical due to the 
smallness of the reduced damping factor a . The whirl frequency ratio for 
synchronous excitation at 1000 rpm is very close to 0.5 due to the smallness of 
the factor a . No significant reduction in damping is found at low frequencies 
due to the low product of recess pressure and compressibility parameter (3. 


The effect of a large compressibility factor 3 on a shallow recess, 
orifice compensated HJB has been considered. The sample bearing is typical of a 
current cryogenic application with liquid hydrogen as the working fluid. The 
geometric and operational parameters for the bearing are given as follows; 

N=8, L=62.58mm(2.463in), D=84.43 mm(3.32in), 

l=b=22 . 776mm(0. 896mm) 

c=0.0254mm( 0.001 in) , Hr/c=5, 10, 15 

d =0. 779mm, 
o 

y=6.279 10" 6 Ns/m 2 , p=51.11 kg/m3, 3=8.86 10" 9 m 2 /N, 

P -P = 20.684 MPa (3000 psig), p =0.5 
s ci ro 

8 = 25000 rpm (416.6 Hz) 

The operating parameters determine the flow regime in the bearing to be 
turbulent, and thus, the results of the present analysis are not valid. However, 
the author (SanAndres , 1989a) has presented a similar analysis for turbulent HJBs 
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with inertia effects. This analysis is more complex and does not provide closed 
form formulae for direct evidence of the compressibility effect on the dynamic 
force coefficients. 

Figures 4 and 5 show the synchronous (u>=Q) stiffness and damping 
coefficients for increasing recess depths and as a function of the 
compressibility parameter 3. For the value of 3 corresponding to liquid 
hydrogen, the calculated break frequency is well above ICrHz and thus, fluid 
compressibility does not affect greatly the direct stiffness coefficients. 

Figure 2 shows the stiffness coefficients to be constant since essentially the 
frequency ratio f is almost zero for the application considered. Figure 5 shows 
the direct damping coefficient to be greatly affected by the compressibility 
parameter 3 and the recess depth. The loss in damping at the operation point for 
LH2 is larger for the largest recess depth considered. 

Figure 6 shows the whirl frequency ratio <p to approach a value equal to 
0.5 as the compressibility parameter tends to zero( incompressible , 3 -»0). As 3 
increases, the whirl ratio <}> also increases denoting a substantial reduction in 
available direct damping. At a value of 3=8.88 10 _ in /N typical for LH2, the 
whirl ratio <j> raises from 0.554 to 0.651 as the recess depth increases from 5 
to 15 times the bearing clearance. 

The results presented clearly show the inconvenience of the proposed HJB 
since the stability characteristics of the bearing are worse than those of a 
plain journal bearing. Further reduction of the recess depth will certainly 
reduce the whirl frequency ratio but will convert the HJB into a plain journal 
hearing and cause a drastic reduction in direct hydrostatic stiffness. Suggested 
fixes that can substantially improve the stability characteristics of the 
bearing are 

a) to reduce the axial extent of the recess so as to reduce the recess volume 
and provide more space for development of the squeeze film damping effect. 
However, shorter recess axial lengths will determine reduced hydrostatic 
stiffness. To reduce the recess circumferential extent is not recommended since 
then the mean flow circumferential velocity will grow rapidly and it may 
determine and effective increase in the whirl ratio. 

b) A roughened bearing surface will reduce substantially the cross coupled 
stiffness coefficient and reduce the whirl ratio <f> veil below 0.5 (SanAndres, 
1989b). Fluid injection opposite to shaft rotation may also reduce the 
cross-coupled effects, although a quantitative analysis is not yet available. 


CONCLUSIONS 

The dynamic behavior of externally pressurized bearings has been shown 
to be highly dependent on the fluid compressibility at the bearing recesses. At 
high frequency excitations beyond a break frequency, the bearing hydrostatic 
stiffness increases sharply and it is accompanied by a rapid decrease in direct 
damping. Also, the potential of pneumatic hammer instability (negative damping) 
at low frequencies seems to be certain in hydrostatic bearing applications 
handling highly compressible fluids. 

The analysis presented has shown that the bearing dynamic force response 
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for the entire frequency spectrum is determined by just 2 parameters, the break 
frequency and the reduced damping factor a . Criteria to avoid stiffness 
hardening and damping reduction by removing tRe break frequency well above the 
desired operating frequency are derived and discussed. At low frequencies, the 
dynamic performance of the bearing is substantially degraded if the reduced 
damping factor is large. Furthermore, the limited stability characteristics of 
hybrid bearings are severely restricted. To avoid this effect it is required for 
the recess volumes(recess depth) to be small. However, for highly compressible 
fluids, improved bearing stability characteristics may demand a reduced limit on 
the maximum external supply pressure to the bearing. Consequently, reduced load 
capacity and stiffness are to be traded off for stability. HJBs with a roughened 
bearing surface could be an appropriate alternative for better bearing stability 
and dynamic performance. 
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TABLE 1 

Dynamic Force Coefficients for Laminar Flow 
Hydrostatic Journal Bearings 

no fluid compressibility effects 
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Figure 1. Geometry of a hydrostatic journal bearing and description of 
coordinate systems. 



Capillar HJB: L/D=,5, 1000rpm, Ps“689500Pa , c=*2.54E-5m 

Figure 2. Effect of excitation frequency on direct stiffness for a 
capillary hydrostatic journal bearing vlth increasing recess 
depths. 
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Direct Damping Coefficient Cxx/Cxxo 


Cxxo-2.5E+6 Ns/m 



Nrecess“8 Frequency ratio, f-W/Wb 

Figure 3. Effect of excitation frequency on direct damping for a capillary 
hydrostatic journal bearing vith increasing recess depths. 

Capillar HJB: L/D— . 5 , 1000rpm, Ps=*689500Pa, c=2.5‘lE-5m 


8-rec H2B 25000rpm 3000 psi, Ine rp - 1 , Iner 1 ■ 1 



Compressibility Parameter B (m2/N) 


ia-Rpm-25000! c=2.54E-5m, Hr/c»5 , 1 0 , 1 5 , 1=22, 7mm 

Figure 4. Effect of fluid compressibility factor 6 and recess depth on 
synchronous stiffness coefficients for a liquid hydrogen, orifice 
compensated, hybrid bearing. 
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S-i'ec H29 25000rpm 3000 psi , 



Compressibility Parameter B (m2/N) 


w-rpm-25000; c=2.S4E-5m, Hrr-c-5 , 1 0 , 1 5 ; 1=22. ?mm 

Figure 5. Effect of fluid compressibility factor 6 and recess depth on the 
synchronous damping coefficients for a liquid hydrogen, orifice 
compensated, hybrid bearing. 
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u»=*Rpm=25000 ; c=2.54E-5m, Hr /c-5 , 1 0 , 1 5 , 1=22. 7mm 

Figure 6. Whirl Frequency ratio + for a LH2 hybrid bearing. Combined effect 
of fluid coapresaibillty and recess depth. 
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